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The parton distribution functions determined by CTEQ at low Q 2 are used as inputs to test the 
D ' validity of the valon model. The valon distributions in a nucleon are first found to be nearly Q 

. independent. The parton distribution in a valon are shown to be consistent with being universal, 

independent of the valon type. The momentum fractions of the partons in the valon add up sepa- 
rately to one. These properties affirm the validity of the valon model. The various distributions are 
parameterized for convenient application of the model. 
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| The parton distributions in proton have been studied extensively in recent years over wide ranges of Q 2 

and x |Q, ||, |j. For example, in the CTEQ global analysis framework |jj the distribution functions have 
been determined by fitting some 1300 data points obtained for various reactions in 16 experiments; over 
30 parameters are used in the input parton distribution functions (PDF's) in the next-to-leading -order 
calculations in perturbative QCD Q. The results on the PDF's at various Q 2 are presented in the form 
Q-i! of graphs available on the web Q. The distribution functions are accurately calculated numerically but 

they are inconvenient to describe analytically. The purpose of this paper is twofold. One is to provide a 
simple parameterization of the parton distribution functions that is accurate to within 5% for 1 < Q 2 < 100 
(GeV/c) 2 and 10~ 2 < x < 1. The other is to provide firm evidence for the validity of a model that can be 
very useful in the study of soft processes in hadronic and nuclear collisions. 

The model under discussion is the valon model fjj. Valons play a role in scattering problems as the 
constituent quarks do in bound-state problems. In the model it is assumed that the valons stand at a level 
in between hadrons and partons and that the structure of a hadron in terms of the valons is independent 
of Q 2 . That is, the property that a nucleon has three valons which carry all the momentum of the nucleon 
does not change with Q 2 . Each valon may be viewed as a parton cluster associated with one and only one 
valence quark, so the flavor quantum numbers of a valon are those of a valence quark. The Q 2 evolution of 
the parton distributions in a nucleon is effected through the evolution of the valon structure, as the higher 
resolution of a probe reveals the parton content of the valons. 

When the valon model was first proposed, the deep inelastic scattering data was not precise enough either 
to rule out the model or to determine accurately the parameters in the model. One may regard the situation 
then as one in which the model satisfies the sufficiency condition for an approximate description of the 
nucleon structure, but not necessary. Now, the experimental data have vastly improved and the PDF's have 
been so precisely determined that the validity of the valon model can be put to a stringent test. That is 
what we intend to do in this paper. Although no model is ever necessary in the mathematical sense, we shall 
show that the concept of valons as constituents of the nucleon is eminently acceptable by virtue of the Q 2 
independence of the valon distribution functions, the universality of the parton distributions in valons, and 
the momentum sum rule of the partons in a valon being satisfied. 

In terms of the parameterization in the valon model the parton distributions are much simpler to describe 
and therefore more convenient to transport to different problems. Our focus will be on soft processes at 
low Q 2 , for which perturbative QCD is inapplicable. The valon picture then provides a systematic way of 
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organizing various contributions to inclusive processes. One such problem is the study of proton- nucleus 
collisions in which the degradation of the momenta of the produced nucleons can be well described in the 
valon model Q . The new parameterizations determined in this paper will affect the details of the model, for 
which we have previously made simplifying assumptions. The application of the results obtained now will 
not be considered here. We mention such applications here only to motivate our concentration in the range 
1< Q 2 < 100 (GeV/c) 2 in this paper. 



II. THE VALON MODEL 



In the valon model we assume that a proton consists of three valons (UUD), which separately contain the 
three valence quarks (uud). Let the exclusive valon distribution function be 

GuuD{yi,V2,ya) = 9 (yilfrTvi 5{vi + y 2 + y 3 - 1), (1) 

where yi are the momentum fractions of the U valons (i — 1,2) and D valon (i — 3). The variable y will 
never refer to rapidity in this paper. The normalization factor g is determined by the requirement that the 
probability for the proton to consist of three and only three valons is one, i.e., 

nl rl—yi rl—yi—V2 

/ dyx I dy 2 / dy 3 GuvD(vi,V2,Va) = 1- (2) 

J o Jo Jo 

Note that the valon distribution function is not defined in the invariant phase space. From Eq. (||) we have 

g = [B(a + l,/3+ l)B(a + 1, a + (3 + 2)]" 1 , (3) 
where B(m,n) is the beta function. The single-valon distributions are 

Gu(y)= [ V dy 2 ( dy 3 G UUD (y,y 2 ,y 3 ) = gB(a + 1,0+ l)y a (l ~ y) a+ ?+\ (4) 

Jo Jo 



Gn(y)= [ J d Vl f " Vl dy 2 G UUD ( yi ,y 2 ,y) = gB(a + 1, a + l)y? (1 - yj 2a+1 . (5) 
Jo Jo 

An essential property of the valon model is that the structure of the proton in terms of the valons is 
independent of the probe. It means that when probed at high Q 2 , whatever the experiment may be, the 
parton distributions in a proton can be expressed as a convolution of the valon distribution and the parton 
distribution in a valon, i.e., 

xu(x,Q 2 )= [ dy[2G u (y)K(x/y,Q 2 ) + G D (y)L u (x/y,Q 2 )], (6) 



x d(x, Q 2 ) = [ dy [G D (y)K(x/y, Q 2 ) + 2Gu(y)L u (x/y, Q 2 )} , (7) 

J X 

where u(x, Q 2 ) and d(x, Q 2 ) are the u and d quark distributions, respectively, and x the momentum fraction 
of the quark. We emphasize that on the right-hand side (RHS) of the above equations the Q 2 dependences 
appear only in parton distributions in the valons, K(z, Q 2 ) and L u (z, Q 2 ), but not in the valon distributions, 
Gu{y) an d Go{y)- We regard this as the defining property of the valon model. In that sense the model is 
analogous to the one in which a deuteron is treated as a bound state of two nucleons; in that treatment a 
high-Q 2 probe resolves the structure of one or the other nucleon without affecting the nucleon wave function 
of the deuteron. 

There are two types of parton distributions in the valons that appear in Eqs. (JsJ) and (^). K(z, Q 2 ) refers 
to the favored partons, i.e., u in U and d in D, whereas L u (z, Q 2 ) refers to the unfavored partons, i.e., u in 
D and d in U. The distribution K(z, Q 2 ) can be further divided into two types 



K(z, Q 2 ) = K NS {z, Q 2 ) +L f {z, Q 2 ), 



(8) 
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where the first term on the RHS is the valence quark distribution (hence, non-singlet), while the second is 
the sea quark distribution for the same flavor type. Since the sea quarks should respect charge conjugation 
invariance, the u and u in the sea (and similarly d and d) have the same distributions, i.e., 



xu(x,Q 2 ) = [ dy[2G u (y)L f (x/y,Q 2 ) + G D (y)L u (x/y,Q 2 )] , 

J X 



(9) 



xd(x,Q 2 ) = dy[G D (y)L f (x/y,Q 2 ) + 2G u (y)L u (x/y 1 Q 2 )]. (10) 

J X 

The valence quark distributions are then 

xu v (x,Q 2 ) = f dy2G u (y)K NS (x/y,Q 2 ), (11) 

(12) 



xd v (x,Q 2 ) = / dyG D (y)K NS (x/y,Q 2 ). 

J X 



In earlier treatment |J-|| Lf and L u have been regarded as identical due to the assumption of the 
symmetric sea. Indeed, putting Lf = L u in Eqs. (||) and (|l"(i| ) results in u(x,Q 2 ) — d(x,Q 2 ), which is a 
necessary consequence of the sea quarks satisfying SU (2) symmetry. However, there is experimental evidence 
H that Gottfried integral J [F^ — F^dx/x is less than 1/3, which is the value expected in the simple quark 
model. Thus we should allow Lf to be different from L u . Indeed, in the valon model we may expect that in 
a U valon the necessary presence of a u valence quark would on the grounds of Fermi statistics make a gluon 
have more difficulty converting virtually into a uu pair than into a dd pair. Hence, Lf should be suppressed 
relative to L u . Since there are two U and one D is a proton, we should expect the overall u to be less than 
d. The data do indicate u < d after integration over x ||. We thus see that the breaking of SU{2) in the 
sea is related to Pauli blocking in the valons. 

III. THE VALON DISTRIBUTIONS 

The valence quark distributions, u v and d v , are given by CTEQ4LQ ||. From Eqs. ([ll]), and ( |l2"| ) we 
see that they are directly related to the valon distributions, Gjj{y) and Go(y), by convolutions with the 
common factor K^s- It is therefore possible to isolate the valon distributions by deconvolution using the 
moments. Let us define 

Gu, D {n)= f dyy^GuMv), (13) 
Jo 

K NS (n,Q 2 )= [ dzz n - 2 K NS {z,Q 2 ), (14) 
Jo 

l 

n-l„ 



u v (n, Q ) = \ dxx n u v (x,Q ), (15) 
Jo 

and similarly for d v in terms of d v . Then by the convolution theorem we have from Eqs. ( pi] ) and (p^), 

u v (n, Q 2 ) = 2Gu(n)K NS (n, Q 2 ), (16) 



d v (n,Q 2 ) = G D (n)K NS (n,Q 2 ). 



(17) 
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It thus follows that 

Gu(n) u v (n,Q 2 ) 



G D {n) 2d v {n,Q 2 ) 



(18) 



If the valon model is valid, then the LHS is Q 2 independent, a property that we can check directly by 
examining the Q 2 dependence of the RHS. Since u v (x, Q 2 ) and d v (x,Q 2 ) can separately be determined from 
H|, we only have to take their moments and calculate their ratio. 

To do the above analysis, we have to set the range of Q 2 to be examined, since the valon model is 
not expected to be accurate for all Q 2 . As we have discussed near the end of Sec. I, the valon model 
has been applied to soft production problems because they involve non-perturbative processes. For hard 
processes at very high Q 2 [Q 2 > 100 (GeV/c) 2 ] perturbative QCD is very successful and there is no need 
to introduce any inaccuracies through the use of a model. We shall therefore limit ourselves to the range 
1 < Q 2 < 100(GeV/c) 2 . This is actually a very wide range for hadronic processes that can involve the 
production of soft particles and semi-hard mini-jets. 

For the range of Q 2 chosen we must use low-Q 2 parameterization of the PDF's. CTEQ4LQ || gives the 
graphs of u, d, s, u, d, and g distributions at any Q evolved from the starting scale at Q 2 — 0.49 (GeV/c) 2 . 
Since u v and d v distributions are not included in the list of PDF's posted, we have to calculate 

u v {x,Q 2 ) = u{x,Q 2 ) -u(x,Q 2 ) (19) 

d v {x,Q 2 ) = d(x,Q 2 )-d{x,Q 2 ) (20) 

from the q and q graphs for the RHS available from the web. We extract the numerical values at up to 
60 points of x values per PDF for three values of Q : 1, 3, 10 GeV/c. From the values of u v and d v thus 
determined, we then compute the moments in accordance to Eq. (|l5|) for n = 2, • • • 9. Denoting the RHS of 
Eq. ( |T8| ) by r(n, Q) we then calculate the ratio of ratios 

R(n,Q)=r(n,Q)/r(n,Q = l) (21) 

relative to Q = 1 GeV/c. Figure 1 shows how R(n, Q) depends on n for Q = 3 and 10 GeV/c. It is evident 
that the dependence is very mild, the maximum deviation from 1 being around 7% at n — 9 and Q = 10 
GeV/c. We regard this approximate constancy of R(n, Q), while Q is increased by an order of magnitude, as 
the first step toward a confirmation of the validity of the valon model. That is, the insensitivity of R(n, Q) 
to Q variation supports our assumption that Gu(n)/GD{n) is Q independent to a degree sufficient for the 
application of the valon model. 

The n dependence of r(n, Q) = u v (n, Q)/2d v (n, Q), as determined from the CTEQ4LQ data, can now be 
used to fix the parameters, a and (3, in the valon distributions. From Eqs. (ji|), (|5|) and (13) we have 



from which follows 



Gu{n) = B(a + n,a + (3 + 2)/B{a + 1, a + j3 + 2) (22) 
G D [n) = B{f3 + n,2a + 2)/B{P + 1, 2a + 2) (23) 



Gu(n) T(a + n)T(p+l) _ 

7(a,/3,n). (24) 



G D (n) T(a + l)T{f3 + n) 
The parameters a and (3 can now be determined by minimizing 

~ 7(q, (3,n) - r{n,Q) 
_7(a,/?,n) +r(n,Q) 



C = J2 



(25) 



where N is the maximum number of moments we extract from the CTEQ4LQ data, which we take to be 
N = 10. Note that n = 1 is excluded in Eq. (|25|), since j(a, j3, 1) = 1 basically due to Eq. (||), and r(l, Q) = 1 
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because there are two u v and one d v . Varying a and (3 in search for a minimum in C results in incredibly 
small C, in the order of 10~ 5 . Wc find 



a = 1.755, /3 = 1.05, 



for Q = 1 GeV/c. 



(26) 



1.545, ,0 = 0.89, 



for Q = 10 GeV/c, 



(27) 



We have not ignored the small Q dependences of a and /3 since the data on r(n, Q) contain some Q depen- 
dence. Moreover, the Q independence of the ratio Gu{n)/GD{n) does not mathematically preclude the Q 
dependences of Gjj(n) and Gn(n) separately. However, the difference between Eqs. (26) and f|27|) is not very 
great, as we can see in Fig. 2, where Gu(y) and Gd{d) are shown [through the use of Eqs. (Q) and (0)] for 
the two extreme Q values, 1 and 10 GeV/c. The difference is insignificant compared to those of the quark 
distributions, one of which is shown in Fig. 3. With the drastic difference between Figs. 2 and 3 in mind, it 
is reasonable to conclude that the essence of the valon model has been verified to the extent that the valon 
distributions exhibit approximate scaling in Q. 

It should be remarked that the values of a and (3 determined above are very different from the ones obtained 
previously. In Ref. |?J the early data of muon [jlOj and electron scattering [ fTl) were used in conjunction with 
a number of theoretical assumptions to yield the values of a — 0.65 and (3 — 0.35. In Ref. Q the modern 
data of CTEQ4LQ were used, but the values a = 0.70 and (3 — 0.25 were obtained due to the assumption 
of a specific form of the sea quark distribution (proven to be grossly inaccurate below). Here in this paper 
we have avoided making any assumption about the sea quark distributions. By extracting the valence quark 
distributions from CTEQ4LQ we have been able to determine Gu{y) and Gd(v) directly. Clearly, the new 
values of a and (3 are more reliable. Further support for their reliability will be given below in connect with 
the quark distributions u and d, for which our previous parameterization in Ref. Q has led to unaccountable 
discrepancies that are unsatisfactory. 

From Eq. (113) one sees that the momentum fractions of the valons are given by the n = 2 moments. From 
Eqs. (p2|), fl23|), (26) and (p7|), one can then calculate the momentum fractions (y), yielding 



(y) u = 0.3644, 
0.3646, 



1 GeV/c, 
lOGeV/c, 



(28) 



(y)i 



At either Q the sum rule 



0.2712 
0.2708, 



Q = 1 GeV/c, 
Q = lOGeV/c. 



2 (y)u + (y) D = i 



(29) 



(30) 



is satisfied identically. We see that the momentum fractions carried by the valons are essentially independent 
of Q and that each U valon carries as much as 1.345 times more than the D valon. 



IV. THE QUARK AND GLUON DISTRIBUTIONS 



Having determined the valon distributions, we can now proceed to the quark and gluon distributions in 
the valons. From Eq. ([!(]) we have 



K NS (n, Q) = u v (n, Q)/2Gu(n, Q) 



(31) 



where wc allow Gu(n,Q) to have its weak Q dependence given by Eqs. ( p6| ) and (27). Using the moments 
u v (n,Q) that we have already calculated from CTEQ4LQ, we obtain the values of Kns(ii,Q) as shown in 
Fig. 4. As expected, Kjys(n,Q) undergoes substantial evolution, especially from Q = 1 to 3 GeV/c. 

To test how good our determination of a and (3 is, we use the R^si™, Q) calculated above in conjunction 
with Go(n,Q) that can be obtained from Eqs. (^|) and (13) so that d v (n,Q) can be computed using Eq. 
(|l7j). Note that this computation of d v (n, Q) requires the knowledge of a and /3, while the computation of 
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d v {n,Q) for the RHS of Eq. @, i.e., r(n,Q), is based on the CTEQ4LQ data for the RHS of Eq. @. 
Our point here is to calculate d(x,Q) from the CTEQ input on u v {n,Q) in Eq. (|3l]), a procedure that is 
made possible by the common factor Kn s(n,Q) in both Eqs. ( |l6| ) and (|l7|). Physically, it means that the 
evolution of quarks in a valon is independent of the flavor of the host valon. 

After d v (n, Q) is obtained by the above procedure that tests the values of a and [3, we then make the inverse 
transform to get the distribution d v (x, Q). This transform can be facilitated by exploiting the orthogonality 
of the Legendre polynomials, the details of which are discussed in Ref. B. Upon the determination of 
d v (x,Q) we can add to it the d(x,Q) distribution from Ref. || and obtain d(x,Q). In Fig. 5 we show the 
u and d quark distributions at Q = 1 GeV/c. The solid lines are the distributions posted by CTEQ4LQ 
The dotted line for xd(x) is what we have computed using the procedure outlined above. Note that its 
agreement with the solid line is excellent. The dotted line for xu(x) is essentially the result from fitting the 
CTEQ4LQ data in the valon model; it merely affirms that the fit is extremely good, so the values of a and 
P are reliable. The result on xd(x) reveals more about the soundness of the model, since it is not obtained 
by fitting, but calculated using the valon distribution Gd{u) and the universality of Kns{ z , Q)- 

For the sea quark and gluon distributions, it is for the convenience of the applications of the valon model 
that we find simple parameterizations of their distributions in a valon. To that end we first write Eqs. (^) 
and (0) in moment form 



u(n, Q) = 2Gu(n)L f (n, Q) + G D (n)L u (n, Q), 



(32) 



d(n, Q) = G D {n)L f (n, Q) + 2Gu(n)L u (n, Q). 
For the strange quark and gluons, we have 

s(n, Q) = \2Gu(n) + G D {n)] L s {n, Q), 



(33) 



(34) 



g(n,Q)= 2Gu(n) + G D (n) L g (n,Q). (35) 



Since u, d, s and g are known from CTEQ4LQ, and Gjj and Gjj known from Eqs. ( p2| ) and (|23|), we can solve 
for Lf, L u , L s and L g . From the result we perform the inverse transform to Lf(z, Q), L u (z, Q), L s (z, Q) and 
L g (z,Q), respectively. These distributions are shown in Figs. 6 and 7, for Q — 1,3, and 10 GeV/c. Being 
in log-log plots, the evolutions due to the changes in Q are substantial, as expected. What we have not 
expected is the drastic difference between Lf(z, Q) and L u (z, Q). For < — ln(l — z) < 0.4, the range of z 
is < z < 0.33. For z > 0.4 we find that Lf(z, Q) < L u (z, Q), at least for Q = 1 and 3 GeV/c. At higher Q 
the evolution can generate more favored sea quarks. Thus at low Q, where sea quarks are few, Pauli blocking 
suppresses the favored sea quarks so much at high z that the unfavored sea quarks dominate. Indeed, our 
calculation of Lf(z,Q) is unreliable for z > 0.4 because from the finite number of moments (n < 10) that 
we have taken the inverse transform generates oscillations in z at high z. There is no such problem with the 
s quark and gluon distributions, as is evident in Fig. 7, since they are not inhibited by Pauli blocking. The 
general property is that all the parton distributions increase significantly at small z, when Q is increased. 

There is no simple way to describe both the z and Q dependences of the parton distributions. In order for 
the valon model to be useful, especially in applications to the study of inclusive cross sections in hadronic 
collisions at low px, an analytic description of each of the parton distributions is needed. For such problems 
only the distributions at Q = 1 GeV/c are relevant, so we fit those distributions by polynomials. The formula 
used for the fitting is 

3 

hxL i (z,Q = l)=^afti, t = -\n{l-z). (36) 

j=o 

The result of the fitting is shown in Fig. 8, for which the values of the coefficients are given in Table I. The 
fits are evidently very good. Thus we have completely specified the PDF's at Q = 1 GeV/c in analytical 
and numerical ways, suitable for transport to problems where such PDF's are needed. 
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TABLE I: Coefficients in Eq. @ 



i 


CO 
Oo 


(0 
«i 


(0 


00 


I 


-2.66 


0.08 


-10.4 


-6.0 


u 


-2.92 


4.0 


-5.95 


-1.4 


s 


-3.30 


-2.4 


2.7 


-1.65 


9 


-0.63 


-3.1 


4.9 


-1.9 



In applications, such as in Ref. m, the moments are more useful than the PDF's themselves. We show the 
moments at Q = 1 GeV/c in Fig. 9, together with their fits. The fitting is done mainly for the convenience 
of applications. The formula used for fitting is 

3 

lnii(n) = -^2bfu j , u = ln(n-l), (37) 

3=0 

where the coefficients are given in Table II. The rapid decrease of the favored quark moments Lf(n) with 

TABLE II: Coefficients in Eq. @ 



i 


b (i) 
°o 


&« 




b (i) 

°3 


f 


4.12 


2.2 


0.2 


0.18 


u 


3.07 


1.5 


0.08 


0.05 


s 


4.21 


1.6 


0.1 


0.02 


9 


0.98 


1.0 


0.05 






increasing n is now very evident, while the other three parton moments have roughly similar n dependences. 

It is important to check the momentum sum rule of the partons. We have seen in Eq. that the valon 
momentum fractions add up to 1; now the parton momentum fractions in each valon must also add up to 1. 
Since Kns(z) and Li{z) are invariant distributions, their moments at n — 2 are their momentum fractions. 
We therefore should have 

K NS (2) + L g {2) + 2 £i(2) = 1 - (38) 

i=f,u.s 

We have J^ats(2) = 0.4707 at Q = 1 GeV/c, and from Table II we can calculate Li{2) = exp [— bj\ ] , yielding 
0.0162, 0.0465, 0.0148, and 0.3754 for i = f,u, s, and g, respectively. According to the LHS of Eq. @ they 
sum up to 1.001, an excellent confirmation of the momentum sum rule. 

As a final item of paremetrizing the moments of the parton distributions, we give here also a formula that 
fits K NS {n, Q) at Q = 1 GeV/c 

3 

\nK NS {n) = -Y,CjU j , (39) 
3=0 

where Cj = 0.753, 0.401, 0.0962, and 0.0555, for j = 0, 1, 2, 3, respectively. The solid line in Fig. 4 shows the 
fit. 

V. CONCLUSION 

We have shown that the PDF's in a proton can be described in two stages: valons in a proton and then 
partons in a valon. The valon distribution functions are essentially independent of Q, while the parton 
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distributions in the valons are Q dependent. The three valons carry all the momentum of the proton, and 
the way that the parton momenta are distributed in a valon is independent of the host valon so long as the 
sea quark flavors are identified as favored or unfavored, instead of by specific flavors like u or d. We have 
found that Pauli blocking significantly suppress the favored quarks compared to the unfavored quarks. At 
Q = 1 GeV/c, the valence quarks carry 47.1% of the proton momentum, while the gluons carry 37.5%. The 
scaling behavior of the valon distributions and the universality of the valon structure give support to the 
valon model as a simple and organized description of the nucleon structure. 

We have determined the parameters in simple formulas that adequately describe the parton distributions 
(and their moments) in a valon at Q = 1 GeV. This is done for the benefit of applications of the valon 
model to low-py hadronic reactions that are not perturbative. Such distributions are needed when the 
parton degrees of freedom are released. If in some nuclear reactions at some energy where the partons do 
not exhibit their dynamical effect beyond the valons, as suggested by Csorgo |L2j for heavy-ion collisions at 
SPS, then the valon distributions are all that is needed. Since the exclusive valon distribution is the absolute 
square of the wave function of the proton in the valon representation, it is also the recombination function 
of valons in forming a proton || . Thus the new values of a and (3 found here affect the calculation of hadron 
production at low px- 

The valon model can, of course, also be applied to other hadrons beside the proton. Although the data 
on the PDF's of the pion or the kaon are not of the same quality as those of the proton, some data on the 
Drell-Yan and prompt photon production initiated by mesons do exist. It will be a natural extension of this 
work to determine the valon distributions in the mesons by using the PDF's obtained from such data. By 
virtue of the universality of the parton distributions in valons, what we have found here from the proton is 
good enough for the mesons also. 

The affirmation of the validity of the valon model makes possible a logical link between the bound-state 
problem of the hadrons in terms of the constituent quarks and the scattering problem in terms of the partons. 
The relationship between the wave functions of the constituent quarks and the valon distributions was studied 
in the context of form factors [|l3|. In view of the new distributions determined here, that problem needs 
to be revisited. On the whole our understanding of the hadron structure problem is enhanced by our study 
here of the modern parton distribution functions in the framework of the valon model. 
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FIG. 1: The ratio of ratios defined in Eq. 




FIG. 3: The u quark distribution functions from Ref. 
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FIG. 5: The u and d quark distribution functions at Q — 1 GeV/c. The solid lines are from CTEQ la]; the dotted 
lines are calculated in the valon model. 
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FIG. 6: The favored and unfavored sea quark distributions in the valon at three values of Q. 
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FIG. 7: The strange quark and gluon distributions in the valon at three values of Q. 
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FIG. 8: Fits of the f,u,s,g distributions at Q = 1 GeV/c, as shown by the dotted lines. 




